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We develop a theory for the vortex unbinding transition in homogeneously disordered supercon¬ 
ducting films. This theory incorporates the effects of quantum, mesoscopic and thermal fluctuations 
stemming from length scales ranging from the superconducting coherence length down to the Fermi 
wavelength. In particular, we extend the renormalization group treatment of the diffusive nonlinear 
sigma model to the superconducting side of the transition. Furthermore, we explore the mesoscopic 
fluctuations of parameters in the Ginzburg-Landau functional. Using the developed theory, we de¬ 
termine the dependence of essential observables (including the vortex unbinding temperature, the 
superconducting density, as well as the temperature-dependent resistivity and thermal conductiv¬ 
ity) on microscopic characteristics such as the disorder-induced scattering rate and bare interaction 
couplings. 

PACS numbers: 74.20.-z, 74.40.-n, 74.62.-c, 74.78.-w 


I. INTRODUCTION 


Disordered superconductors, the superconductor- 
insulator quantum phase transition, and quantum trans¬ 
port through superconducting nanodevices remain a field 
of intense research over the past 50 years. On the 
one hand, superconducting electronic components are 
of great technological interest in connection with a va¬ 
riety of possible applications. On the other hand, 
the interplay of two most outstanding manifestations 
of quantum coherence—superconductivity and Anderson 
localization—determines the attention attracted by these 
systems in the context of fundamental physics research, 
see, in particular, Refs. 1 3 and references therein. 

Superconductivity, i.e., the phenomenon of frictionless 
transport and perfect diamagnetism, is a consequence of 
long-ranged correlations of the complex order parameter 
A (a:) in a theory of charged particles: 


(A*(x)A(0)> x 


e X ^ A , normal state; 
|(A)| 2 , superconductor. 


( 1 ) 


Here, is the correlation length. In two spatial dimen¬ 
sions at finite temperature true long-range order is not 
possible in view of the Mermin-Wagner theorem. In this 
case, one resorts to the following weaker definition: 


homogeneous. (ii) Strong phase fluctuations of <f> = 
arg((A(a:)}) are suppressed due to sufficiently large 
phase rigidity. 

As a consequence of these conditions, two different 
mechanisms driving the transition between the supercon¬ 
ducting and the normal state are often distinguished: (i) 
the expectation value (A(x)) vanishes across the transi¬ 
tion; (ii) the expectation value (A(x)) + 0 is locally finite, 
but the phase rigidity vanishes across the transition. 

The first of these two mechanisms is sometimes referred 
to as “fermionic” scenario. It includes the Bardeen- 
Cooper-Schrieffer (BCS) theory 4 and related theories. In 
contrast, in the second “bosonic” mechanism, the phase 
fluctuations of preformed Cooper pairs drive the transi¬ 
tion; typically the fermionic spectrum displays a pseudo¬ 
gap even in the normal state. 

A particularly important representative of bosonic the¬ 
ories is the Berezinskii-Kosterlitz-Thouless (BKT) tran¬ 
sition 5,6 occuring in two dimensional (2D) films. In a 
system with broken U(l)-symmetry the phase degree of 
freedom corresponds to a Goldstone boson, the latter 
being described by the following U(l) nonlinear sigma 
model (NLcrM) action: 

Su(i) = ~ f (V^) 2 ■ (3) 

7T J x 


(A*{x)A m 


e X ^ A , normal state; 
l/x* 1 , superconductor. 


( 2 ) 


In the last equation, the power 77 takes values 0 < 77 < 1. 

Typically, the following two sufficient conditions are 
fulfilled in a superconductor: (i) The modulus of the 
expectation value |(A(a:))| is nonvanishing and nearly 


Here and below we use a short-hand notation f x for the 
spatial integral over the 2D system / d 2 x. In a clean 
superconductor the phase rigidity is given by K/ir - 
n s /2m*T, where the density of Cooper pairs of mass 
in* = 2m is denoted by n s . We set both the Boltzmann 
and Planck constants to unity, ks = ft = 1, throughout 
the paper. 
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Since the manifold of Goldstone bosons is flat, the 
U(l) NLctM, Eq. (3), is not renormalized perturbatively. 
In other words, the theory is Gaussian on the pertur¬ 
bative level. However, the inclusion of nonperturbative 
effects (vortices) leads to a two-parameter renormaliza¬ 
tion group (RG ). 7 If the stiffness is less than the critical 
value It* = 1, vortices proliferate and the system can not 
sustain phase coherence. Contrary, for K > A'* there is 
a regime where vortices are tightly bound into dipoles 
and the system is superconducting. The critical stiffness 
defines the BKT transition temperature Tbkt- If we as- 
surne n s to be determined by the BCS expression, Tbkt 
is parametrically close to the BCS transition temperature 
Tbcs : 


The goal of this paper is to develop a theory of 
the BKT transition in disordered superconductors that 
takes into account the localization and the Coulomb- 
interaction physics. Let us first remind the reader about 
key previous works that, while having strongly advanced 
the understanding of the disordered BKT transitions, did 
not include effects associated with the localization and 
the Coulomb interaction. 

In full analogy to the clean case, the disordered BKT 
transition is determined by A'* = 1, with the only differ¬ 
ence being 


AT = 


ngp A 
16 T 


, A 
tann — 
2 T 


( 8 ) 


Tbkt = Tbcs (1-4 Gi). (4) 

Here, Gi = Tbcs/i* « 1 is the Ginzburg-Levanyuk num¬ 
ber 8 for clean 2D superconductors 9 and /i denotes the 
chemical potential. 

In this paper we concentrate on disordered systems, 
where the elastic scattering rate 1 /r of the electrons sat¬ 
isfies the condition 


Tbkt « 1/t « p . (5) 

According to “Anderson’s theorem”, 10,11 the critical tem¬ 
perature T c with and without disorder are equal T c = 
Tbcs in the absence of electron-electron interaction in 
singlet or triplet channels. The “Anderson theorem” 
discards, however, two important quantum-interference 
phenomena that can dramatically affect the supercon¬ 
ductivity. The first one is the disorder-induced Anderson 
localization . 12 The localization is a natural antagonist 
of superconductivity, and a competition between them 
leads to emergence of superconductor-insulator transition 
in 2D systems. Second, the interelectronic Coulomb in¬ 
teraction can drastically change the critical temperature 
T c , see Ref. 13 for an early work in this direction. The 
analysis of RG equations for the interacting NLcrM of dif¬ 
fusive soft modes 14-18 predicts suppression of T c , which 
is governed by the following expression 19-22 



In this formula 7 *°^ = , z *, 1 -r. The result is valid in 

the limit of sufficiently large dimensionless Drude con¬ 
ductance go = 2/(7 ito) » 1 in units of e 2 //i. Finally, 
for short range interaction, the mean-field transition 
temperature is predicted to increase in the presence of 
disorder , 3,23-25 (provided | 7 c°^ | « t) 



The physical mechanism behind this phenomenon is wave 
function multifractality, which leads to an enhancement 
of the matrix elements of interaction . 2 ' -28 


in the case of disordered films. As a consequence, in the 
disordered case , 26 Tbkt is again given by Eq. (4) with 
Gi now being 9 


Gi = 


7C(3) 

ir 3 g D 


« 1 . 


(9) 


Most prominently, the superconducting transition 
manifests itself in the temperature dependence of con¬ 
ductance. There is a parametrically small tempera¬ 
ture window close to Tbkt within which the normal 
state resistance rapidly drops to exponentially small val¬ 
ues. Above T c , fluctuating Cooper pairs lead to an 
enhancement of conductivity, via Aslamazov-Larkin, 
Maki-Thompson 30,31 and density of states 32 (DOS) 
corrections . 9 In the superconducting state, the resistance 
is determined by vortex excitations in the order parame¬ 
ter field . 33 Due to the Joseplison relation <j> - 2eV a finite 
resistance (steady voltage V) requires a phase relaxation 
mechanism. In the 2D case close to the thermodynamic 
transition, the phase relaxation rate is determined by 
the vortices traveling across the system perpendicularly 
to the current direction . 34,35 Above the vortex unbind¬ 
ing temperature, the latter proliferate and the following 
crossover formula for the resistance p near T c was sug¬ 
gested by Halperin and Nelson in Ref. 36 


P(T) = t n 


1 


1 + 



€hn (T) \ 
£(7bkt) / 


( 10 ) 


The numerical constant A v is fixed to A v = tt 4 /56C(3) * 
1.45, so that Eq.(10) interpolates between the vortex¬ 
generated resistance in the superconducting state and 
the Aslamazov-Larkin contribution above T c . Here the 
length £hn(T) is given by 


JaL) ( 11 *) 

0 \ V tbkt ) 


and provides an interpolation between the BKT coher¬ 
ence length 

£hn £bkt(T) = ^ - T) ex P(^VGT/ tbkt) (Hb) 
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in the limit tbkt « Gi and the mean field behavior 

£hn -»■ £gl(T) = -y===, (He) 

VTbkt 

in the opposite limit, tbkt » Gi. The quantity b is a fit¬ 
ting parameter of order unity. 37 Further, £(Tbkt) is the 
Ginzburg-Landau coherence length given by £(Tbkt) = 
£(0 )/\/Gi, and tbkt denotes a relative distance to the 
BKT transition, 

tbkt = (T - T bkt )/T bkt . (12) 

The observation of the BKT transition in supercon¬ 
ducting films is a longstanding experimental challenge. 
The major difficulties are due to (a) the parametrically 
small window T c - Tbkt « T c and (b) finite size effects. 
Both the finite thickness 38 and the finite area 3 '’ 39 of the 
film lead to an infrared cut-off in the logarithmic vortex- 
antivortex interaction and may thus obscure the BKT 
physics. Defining features for the experimental observa¬ 
tion of the BKT transition are (i) a universal crossover 
function In (p(T)/to) <x 1/x/tbkt, cf. Eq. (10); (ii) a 
jump from linear to nonlinear resistance (V oc /“) with 
a = 3 right at Tbkt; (iii) a crossover from sublinear to su- 
perlinear magnetoresistance. For early work in this area, 
we refer the reader to Refs. 26,40-42 and to the review 
43. Recent years have witnessed a strong rise of exper¬ 
imental interest to this phenomenon, with observations 
of the BKT transition in superconducting films made of 
various materials, including indium oxide, 44-46 MoGe, 46 
titanium nitride, 4 ' niobium selenide, 48 lead, 49 niobium 
nitride, 0-52 and iron selenide. 53 

It is important to emphasize that the experimen¬ 
tally determined Tbkt may differ very substantially from 
Eqs. (4) with an input from Eq. (9); see, in particu¬ 
lar, Ref. 47. The reasons for the insufficiency of the 
theory leading to Eqs. (4) and (9) is closely related to 
those of the failure of Anderson’s theorem, see a discus¬ 
sion and references above. These equations are based on 
the expression (8) for the superconducting density that 
does neither take into account the strong renormaliza¬ 
tion of conductivity due to interference and interaction 
effects, nor the renormalization of T c , see, in particular, 
Refs. 54,55. As a consequence, Eqs. (4) and (9) are not 
sufficient to correctly predict Tbkt as a function of mi¬ 
croscopic parameters encoded in r oc g D and Tbcs- 

We are thus facing the following questions: (i) What is 
the superconducting stiffness including the disorder- and 
interaction-induced corrections? (ii) What is the vortex 
unbinding temperature for homogeneously disordered su¬ 
perconducting films? (iii) What is the temperature de¬ 
pendence of resistivity? The goal of the present paper 
is to develop a theory that answers these questions. As 
we will show, this requires an implementation of a strat¬ 
egy that allows one to go from the Fermi-liquid theory 
at relatively high energies to the low-energy U(l) theory 
through a sequence of intermediate-scale effective field 
theories. We will assume throughout the paper, that 


the dimensionless resistance is small for all temperatures, 

p{T) « 1. 

The paper is structured as follows. Section II is de¬ 
voted to the determination of the vortex unbinding tem¬ 
perature Tbkt- As will be explained in detail in Sec. II A, 
our formalism is based on the consecutive use of Fermi 
Liquid (Sec. IIB), diffusive NLctM (Sec. IIC), Ginzburg- 
Landau (Sec. IID) and U(l) NLctM (Sec. HE) theories. 
The subsequent Section III is devoted to the temperature 
dependence of resistivity of the metallic film close to the 
superconducting transition. We conclude the paper with 
a summary and outlook. The most technical details of 
our calculations are delegated to a number of appendices. 

II. FIELD THEORY OF DISORDERED 
SUPERCONDUCTORS 

This section is devoted to the theoretical framework 
of the present paper and of disordered superconductors 
in general. It is instructive first to get a feeling for 
the energy and respective length scales in the problem 
(Sec. II A). From the hierarchy of length scales, the 
Fermi liquid (Sec. IIB), the diffusive NLctM (Sec. IIC), 
the Ginzburg-Landau theory (Sec. IID) and the U(l) 
NLctM (Sec. HE) appear as a sequence of theories. This 
section is structured following this hierarchy and in each 
subsection IIB-IIE we discuss the derivation and the per¬ 
turbative renormalization of the corresponding theory. 

The physical observables discussed in this section are 
the vortex-unbinding temperature Tbkt in a disordered 
superconductor as well as the superconducting density, 
see Sec. HE. 


A. Strategy and hierarchy of length scales 

We will be interested in temperatures close to the vor¬ 
tex unbinding transition: |tbkt| « 1. Our calculations 
are controlled in the limit, when the normal state con¬ 
ductance close to the transition is large: g » 1. In this 
case, Tbkt turns out to be parametrically close to the 
mean field transition temperature Tmf- As we will dis¬ 
cuss below [see specifically Sec. IID 4 and Eq. (64)], the 
latter is close but remains below the critical temperature 
T c associated with the BCS-like instability. Note that in 
general, g + go and T c + Tbcs, see, e.g., Eqs. (6) and (7). 

On length scales larger than the elastic mean free 
path l = vft (vf is the Fermi velocity) we will self- 
consistently associate a length scale Le to an energy E 
by Le = \/D(Le)/E and reversely define El = D(L)/L 2 . 
Here, D(Le ) is the diffusion constant at the scale Le and 
D(l) = Vpl/2 ■ Thus, the three length scales L T ~ L Tbkt ~ 
T'T’mf are close to each other in the regime under consid¬ 
eration. This regime further implies a hierarchy of length 
scales that we shall expose in this section. This hierar¬ 
chy is associated with a step-by-step quantum to classical 
crossover governed by subsequent freezing of excitations. 
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Fermionic Landau quasiparticles with well-defined mo¬ 
mentum are good excitations only on length scales 
shorter than the mean free path l, see Sec. IIB. On longer 
length scales they are “confined” in diffusive soft modes: 
noninteracting diffusons and cooperons as well as inter¬ 
acting bosonic modes in the Cooper and particle-hole in¬ 
teraction channels. In general, their interplay leads to 
strong renormalization of the conductivity and of the 
mean field transition temperature Tmf- In our approach, 
this effect will be captured by the RG technique applied 
to the interacting, diffusive NLctM. The RG stops at a 
length scale Lp a parametrically close to Tt mf (more de¬ 
tailed explaination can be found in Sec. IIC). The only 
modes that remain at larger distances are static fluctu¬ 
ations of the order parameter field A(x). Thus, at the 
scale Lt c we derive the Ginzburg-Landau (GL) free en¬ 
ergy functional weighting these fluctuations, see Sec. IID. 
Our derivation, which contains terms beyond leading or¬ 
der in 1/g, also yields random fluctuations of the coeffi- 
cents in the GL functional. For T < Tmf the Higgs field 
[fluctuations in the modulus of A(x)] is gapped on the 
scale of the GL coherence length £ = £gl(T) » Tt mf - In 
the symmetry broken state, only phase fluctuations of the 
order parameter field are important on scales exceeding 
£. Thus, at the scale of the coherence length we derive 
the disordered model of phase fluctuations, Sec. HE (see 
also Ref. [56] for a related study). Disorder terms turn 
out to be RG-irrelevant and can be taken into account in 
a perturbative manner. This eventually leads to an ac¬ 
tion analogous to Eq. (3), but with a renormalized stiff¬ 
ness K. For Tbkt < T < Tmf, the theory also predicts 
a renormalized coherence length £bkt > £ beyond which 
phase correlations decay exponentially. 

All in all, we find the following hierarchy of length 
scales (the Fermi wavelength is denoted by X F ): 

X F < l < L t < £ < £bkt- (13) 

At each intermediate length scale, a certain “micro¬ 
scopic” theory ceases to be the appropriate description 
and we derive an emergent effective theory valid at longer 
length scales: 


scale 

“microtheory” 

emergent theory 

l 

Fermi liquid - 

* diffusive NLctM, 

Lt c 

diffusive NLctM - 

* GL theory, 

s 

GL theory - 

* U(l) NLctM. 


A summary of length scales, relevant excitations and ef¬ 
fective theories applicable to the various regimes is given 
in Fig. 1. 

B. Fermi liquid theory 

The main statement of Landau’s Fermi liquid the¬ 
ory is that, in the absence of spontaneous symmetry 
breaking, the low-energy excitations of a strongly corre¬ 
lated fermionic system are fermions (Landau quasiparti¬ 


cles) with the same quantum numbers as the free par¬ 
ticles. Their decay rate is small as compared to the 
Fermi energy. In field-theoretical language, this state¬ 
ment means that the exact electronic Green’s function 
(i.e. two point correlator) can be shown to contain a sin¬ 
gular part (quasiparticle pole) with a weight 0 < a < 1 
and an additional regular contribution. 57 

Another particularly important quantity in the the¬ 
ory of strongly interacting fermions are the four point 
correlators. 58-60 These implicitly define the full interac¬ 
tion amplitudes. The latter are subdivided into different 
channels of small energy-momentum transfer according 
to their tensor structure in spin space. For the problem of 
disordered, interacting, spinful fermions we concentrate 
on the static part of the Cooper singlet (T c ), particle- 
hole singlet (T s ) and particle-hole triplet (T t ) scattering 
amplitudes and keep only their zeroth angular harmonic 
(s-wave). 16,18 In this paper, the quasiparticle residue a is 
absorbed into the definition of fermionic fields and scat¬ 
tering amplitude. 

Even though our goal is to describe a superconduct¬ 
ing system close to and below Tmf, be. in the symmetry 
broken phase, it is appropriate and justified to describe 
it using the Fermi liquid theory at the smallest length 
scales L < l. The fermionic excitations at these scales 
do not “know” about the fact that at larger length scales 
they will eventually form coherent Cooper pairs. In other 
words, if the system had a linear dimension L n < l, the 
associated size-quantization energy scale would by far ex¬ 
ceed the superconducting gap. 


C. The diffusive NLctM 

In the previous section, we explained that the Fermi- 
liquid description is appropriate for disordered supercon¬ 
ducting films at smallest length scales, i.e., those less 
than the elastic mean free path l. In our hierarchy of 
length scales we now reach the next level characterized by 
scales exceeding l. The effective field theory that emerges 
in this regime is the NLctM of diffusive interacting soft 
modes (for review see Refs. [16,17]). 


1. Normal state NL<jM 

Upon inclusion of sufficiently weak disorder (in the 
sense go » 1) static quantities, such as the static in¬ 
teraction amplitudes, remain unchanged even at scales 
L » l. They are determined by scales much shorter than 
the mean free path. On the contrary, the dynamical prop¬ 
erties of the system are altered, as the retarded-advanced 
ladders consisting of two Green’s functions acquire a dif¬ 
fusive pole. 16,18 It is possible to describe the diffusive dy¬ 
namics by means of the interacting, diffusive NLctM. 14,15 
In the normal state, the path integral representation of 
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FIG. 1: Hierarchy of length scales for T - Tbkt < Tmf- The scales Lt x and Lt c are introduced for technical reasons, see 
Sec. IIC, below. Both of them and Tt bkt are parametrically close to Lt mf ■ In the category ‘excitations’, a colored bar indicates 
the regime of importance of the various excitation modes. Analogously, in the category ‘theories’, a colored bar indicates the 
regime, where a certain theory provides an appropriate description. A general explanation of the hierarchy of length scales can 
be found in the main text, Sec. II A. 


the partition sum 

Z = f VQexp(-S[Q]) (14) 

is governed by the following action: 

S = S <r + S<£ + Sf£ + sg>, (15a) 

with 

S a = -^ f tr[(VQ) 2 ] fti[ZQ], (15b) 

O x x 

f tr [I^troQ] tr [I- n t r oQ ], (15c) 

4 a,n Jx 

r= 0,3 

s£? = ^r t £ f tr [I%t rj Q] tr [r n t rj Q] , (15d) 

4 a,n Jx 

r= 0,3 
3= 1,2,3 

= ^r c Y, [ tr [KUoQ] tr [Kt r0 Q ]. (15e) 

4 a,n Jx 

r= 1,2 

Here, the symbol ‘tr’ denotes summation over all inter¬ 
nal matrix indices. We are interested in systems where 
time-reversal symmetry is fulfilled. If the system is addi¬ 
tionally spin-rotational invariant (which corresponds to 
class AI in the classification of non-interacting systems), 
the Q matrices are symplectic, traceless, and have non¬ 
trivial structure in replica, Matsubara, spin and Nambu 
spaces: 

Q = = Q _1 = t 12 Q T t r 2 , tr Q = 0. (16a) 

In the absence of spin rotational symmetry (noninteract¬ 
ing class All), the Q matrices are orthogonal, 61 trace¬ 
less, have nontrivial structure in replica, Matsubara, and 


Nambu spaces, and are proportional to the identity ma¬ 
trix in spin space: 

Q = = Q~ x = t w Q T t w ex l a , tr Q = 0. (16b) 

We use the convention t r j = r r <S> aj where r r = (l T ,f) 
are the identity and the Pauli matrices in Nambu space, 
while oj = (l CT ,tf) are those in spin space. Here and 
throughout the paper we use a convention in which a, /3 = 
1,..., Nr denote replicas and m,n = ~N ' M ,..., N' M - 1 
Matsubara indices associated to fermionic frequencies 
e n = 7rT(2?r + l). The following matrices, which are trivial 
in Nambu and spin spaces, have been introduced: 62 


A “ 0 

nm 

= sgn (n) 5 a ^S nm , 

(17a) 

?“/3 
c nm 

= e n S a0 S nm , 

(17b) 

(Ta 0 \ a P 

V n o / nm 

- a^o^A^o^A 

— u u u n-m,n o > 

(17c) 

<0. £ 
s g 

o o 
8 g 

= S a ° a S ao0 6 n+m+hno . 

(17d) 


The coupling constants of the NLuM are the dimen¬ 
sionless conductivity g (bare value go), the static inter¬ 
action amplitudes Tj (i = s,f,c), and the prefactor Z u 
of the frequency term, which is related to the renormal¬ 
ization of specific heat. Note that the latter does not 
flow in the noninteracting case and keeps the bare value 
Z = nis/A. In the presence of long range Coulomb in¬ 
teraction the NLcrM is “JF-invariant”. 63 Essentially, this 
means electrostatic gauge invariance (i.e. invariance un¬ 
der time dependent but space independent phase rota¬ 
tions) and fixes Z u + T s = 0. In the present convention, 
attraction in the Cooper channel implies T c < 0. 

For the sake of a better readability, we omit gauge po¬ 
tentials in Eq. (15) and in the rest of the paper (except 
App. B). Thus, we formally treat a neutral superfluid. 
As we explain in Sec. Ill C, in the truly 2D limit all pre¬ 
sented results hold equally for a charged superconductor. 
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2. NLoM in the superconducting state 

While in the normal state diffusive fluctuations are as¬ 
sociated to smooth variations around the noninteracting 
saddle point solution Q - A, in the symmetry broken 
phase the true saddle point is a function of the supercon¬ 
ducting gap. This result can be directly derived from the 
microscopic theory, see Refs. 64-67. However, in such a 
derivation the disorder is taken into account at the level 
of Anderson’s theorem only. Here we derive the NLctM 
for the symmetry broken phase directly from the interact¬ 
ing normal state NLctM, Eqs. (15). Since this procedure 
can be accomplished at any scale, this allows us to go 
beyond Anderson’s theorem. 

Upon Hubbard-Stratonovich decoupling in the Cooper 
channel the action becomes 

ShsB.A] - ^Tr[(VQ) 2 ]-2Z„Trt3 , + ,( “ q A 1 

+ (18) 

n Jc- 1 - an u x 

The symbol ‘Tr’ includes both the trace operation in 
matrix space and spatial integration. The complex 
Hubbard-Stratonovich field A is defined by 

A=^A a>m L“. (19) 

a,m 


Then A has the following structure in Matsubara and 
Nambu spaces: 



In the limit A -> 0 the usual diffusive form is restored: 
A ^ A = diag(l,-l). 

Returning to the generally complex A Q , it is possible 
to perform a Bogoliubov transformation, 64 

Q-UiqU A , (24) 

such that the saddle point of the rotated field is again 
q = A. In the parametriztion A Q = \A a \e l ^ a the unitary 
rotation matrix is 

Ua= y pa e i^T z ( cosC T v sinC \ 

n>0,a W \-T y sin ?/>“ cosV>“ / (A) 

(25) 

We introduced the energy dependent rotation angle 


We also introduced the notation = T c /Z u . Since we 
expect static s-wave superconductivity, the Hubbard- 
Stratonovich field can be assumed to take the form 
A ai „ = A a S n fi on mean field level. We will refer to A Q 
as “order parameter field”. 

Variation of Eq. (18) with respect to A* 0 leads to the 
gap equation 






( 20 ) 


Using the constraint (16a) [or (16b)] on the Q-matrix, 
we find that the second term in (18) modifies the saddle 
point Q = A to Q = A with 


a = E P\° 

n>0,a: 


|e ra ]A z 

4 + |A 0 


By PA we denote a projector on replica a and on a block 
with given modulus of the Matsubara frequency, matrices 
A x> y jZ are Pauli matrices in this space. 

On the saddle point level, it is possible to choose the 
order parameter field real and equal in all replicas 



After integration of fluctuations 5A a>n = A Qj „-A Q (5 ni o 
around the mean field solution we arrive at the NLctM 
describing the system at T < Tmf- It has the standard 
structure exposed in Eq. (15a). The terms = 

S&WaQU*] and S£>[Q] = S^[UlqU A ] are deter¬ 
mined by Eqs. (15c) and (15d) respectively. The Cooper 
channel interaction term, Eq. (15e) is slightly modified, 
as the static mean field solution is subtracted from Q: 

S int = VT r c E [ tr[Kt r0 Q M ]tr[LZt r0 Q M ]. (27) 

4 a,n Jx 
r= 1,2 

with = Q - AS n (j. The major modification concerns 
the dynamic part of the action, i.e. Eq. (15b), which is 
most conveniently written in the rotated basis, 

S *M = J^Tr [(Vg) 2 ] - 2Z u Tr[qi\. (28) 
We have introduced the matrix 
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3. Saddle point equation 


4- Renormalization group flow 


Using Eq. (23), the gap equation (20) becomes 


A 

T 


oo 

- 7 c ttA Y, 

n =-oo 


l 



r x tanh ; ^ 

/ de " =• 
Ja sjffl _ A 2 


(30) 


This equation has the structure of the standard BCS 
equation. The symbol “=” indicates equality of sum and 
integral upon appropriate ultraviolet (UV) regularization 
at the scale A. The solution of this equation determines 
a transition temperature 


In the normal state, the noninteracting diffuson of a 
quasiparticle pair with Matsubara frequencies e ni , e „ 2 
(ni > 0 , ri 2 < 0 ) and replica indices a,0 is determined 
by the following diffusion propagator (see App. A) 


Me)] 


a(3 

ni,ri2 


D 

D q~ + Cni — C n2 


(36) 


The diffusion coefficient D is determined by the coupling 
constants of the NLcrM via D = g/WZ^. In the super¬ 
conducting state this propagator becomes 


T c = Aexp(l/ 7 C ). (31) 


All standard implications (e.g. A (T), the DOS etc.) im¬ 
mediately follow analogously to the BCS case. How¬ 
ever, it is important to keep in mind, that in general 
the Cooper channel interaction parameter y c = 7 C (A) is 
strongly scale dependent and influenced by disorder and 
interactions in the particle-hole channels. 

In the simplest approximation, the NLcrM is analyzed 
at bare level. The diffusive UV cut-off is A ~ 1/r and y c is 
replaced by 7 C (Z). It already includes the ballistic renor¬ 
malization from scales between the Debye wavelength A jj 
and the mean free path l. Within this simplified analysis, 
T c is determined by the energy scale of divergence of the 
solution 


(<?)] 


a.j5 

ni,ri2 


D 


Dq 2 + e 


a 
n i 


+ £ r 


(37) 


However, in the interval of length scales l < L < Lp « La 
the effect of the superconducting gap on Eq. (37) is neg¬ 
ligible. Therefore, in this interval, it is legitimate to con¬ 
struct a perturbative RG around the normal state sad¬ 
dle point A. (Fluctuations are too fast to resolve the 
difference between A and A.) Thus, in this regime of 
length scales, the renormalization of parameters in su¬ 
perconducting state is dictated by the RG equations of 
the normal state NLcrM, Eq. (15). These have the fol¬ 
lowing form 


7c ^ [7c(0] 1 - In Epr 

of the BCS like 68 RG equation 


djc 

dy 


- 2 7 


2 

C 5 


(32) 


(33) 


where 


dt 

~r = Pt (U 7s; 7t) 7c) , 
dy 

(38a) 

~T = Pis U>7s>7i,7c) , 
dy 

(38b) 

d lt 0 x 

~r = Pit (C 7s; 7t) 7c), 
dy 

(38c) 

d 7c 0 n \ 

— = p 7c {t, 7 s ,7t,7c) • 
dy 

(38d) 


y = -ln(E L T)/2 (34) 

is the logarithm of the running RG scale. The transition 
temperature appearing in Eq. (31) simply becomes the 
BCS transition temperature T c = Tecs with 

1 1 1 1 
T c = — eTTto = WD e 1 ' (A D) = T c e Tc(iT c ) , (35) 

r 

where u>d is the Debye frequency. This equation is a 
restatement of Anderson’s theorem. For BCS supercon¬ 
ductors, 7 c (Ab>) determines the bare phonon mediated 
interaction vertex. 

In general, disorder and interaction effects modify the 
RG equation (33). Then, the transition temperature T c 
(defined via the scale where y c diverges) may strongly 
differ from Tbcs- 


It is worth stressing that, as a consequence of dimen¬ 
sional analysis, the RG equations can be written in terms 
of reduced coupling constants 7 j = Ti/Z u (with i = s, t , c). 
The particle-number conservation implies that the com¬ 
bination Z^ + T s is not renormalized. Therefore 

d\n(Z w ) (t,7s,7t>7c) {no \ 

— 1 — = - n -■ 

In Ref. 15, the beta functions, Eqs. (38), were derived 
to the lowest order in t and j c . Recently , 3 three of us ex¬ 
tended the results of Ref. 15 by deriving one-loop (lowest 
order with respect to t) beta functions which are for¬ 
mally exact in the Cooper channel interaction constant 
7 C . Here, we only quote the result and refer to Ref. 3 for 



more details: 


Pt = t 2 [l + /(7s) +3/(7t)-7c], (39a) 

Pis = (! + 7s) (is + 3 7t + 2 7c + 4 7 ^) , (39b) 

Pit = ~\( l + 7t) [7s -lt~ 2 7 c (1 + 2 7* - 27 c )] , (39c) 
Pic = -^[ (! + 7c) (7s - 3 7t ) - 27 1 + 4y 3 
+ 6 7 c (7t - In (1 + 7t)) ] - 2 7 c • (39d) 

The function fix) entering Eq. (39a) is given by 

f(x) = 1--—— ln(l + x). (40) 

x 


Equations (39) are appropriate for a system with spin- 
rotation invariance. The first term “1” in the square 
bracket of Eq. (39a) describes the weak-localization ef¬ 
fect, which originates from disorder and is unrelated to 
interactions. The last term “- 27 c” i n Eq. (39d) rep¬ 
resents the Cooper instability, which is also present in 
clean systems . 68 All other terms stem from the interplay 
of disorder and interactions. 

In the case of a system with strong spin-orbit coupling 
the following modifications to Eqs. (39) occur. First, 
one should replace the weak-localization by the weak- 
antilocalization effect, i.e. “ 1 ” in the square bracket 
of Eq. (39a) by “-1/2”. Second, the triplet channel is 
gapped out, so that Eq. (39c) should be discarded and 
terms containing should be removed from the remain¬ 
ing equations. 


5. Range of applicability of perturbative RG 

As we have already stated, per definition 
Lt-*Lt c 

7 c {Lt) —* - 00 . Thus, close to T c , one may 
be tempted to keep only the leading powers of 7 C in 
Eqs. (39). However, one should keep in mind that the 
RG Eqs. (39) were derived in the one loop approximation 
(i.e. perturbatively in resistance t). An inspection of 
the perturbative series in the vicinity of T c (where | 7 C | 
is large) indicates 3 that the actual parameter of the 
expansion in this region is t| 7 c |. Thus, close to T c , 
the RG equations are only applicable for energy scales 
El £ Tx , where Tx is defined by 


\lc{L Tx )t{LTx)\ = 1. (41) 


Therefore, disorder-induced corrections are subleading 
with respect to the dominant Cooper-instability term 
within the range of applicability of Eqs. (39). There¬ 
fore, close to T c , the Cooper channel coupling constant 
diverges as 


7 c{Lt) 


L r x'—*L 7 


T r 


( Lt Z Lt x < Lt c )• (42) 


We remind the reader that T c (and thus Lt c ) as given 
by the RG equations (39) is in general strongly renor¬ 
malized as compared to T B cs> see Eqs. ( 6 ) and (7) of the 
introduction. 

We expect that the full RG equations (38) contain a 
line of attractive fixed points characterized by 7 c (Lt 0 ) = 
-00 but finite Z u = Z u {Lt 0 ) and t = t(LT c ) < tsiT ~ 1- 
This expectation is based on the fact that Z u and t deter¬ 
mine the coefficients of the GL functional (see Sec. IID 
below). The line ends in the point with coupling con¬ 
stant igix which represents the quantum critical point 
of the superconductor-insulator transition. While it is 
worth emphasizing the difference between the physical re¬ 
sistance p and the NLcrM coupling constant t (see Sec. Ill 
below and Ref. 3) at this quantum critical point we ex¬ 
pect p = t = tsrr- 

We now present an estimate for the difference between 
T c , t(L Tc ), Z u (L Tc ) and T x , t(L Tx ), Z u {L Tx ), respec¬ 
tively. We will assume that (a) there is a region of fixed 
points, (b) 7 “ 1 -a c ln(T c /T) as Lt -*• Lp c (a c ~ 1 ), (c) the 
crossover of RG flow from Eqs. (39) to the fixed point re¬ 
gion is encoded in a power series of the parameter | 7 c t| 
entering the beta functions, Eq. (38). Mathematically, 
these assumptions mean 


Pt = f/t(l7cf|), 


(43a) 

Pis = - 2 t 1 (l + 7 s )/ 7 s (l7cf|), 

(43b) 

Pic = - 2 lcfic(hct\), 

(43c) 

with interpolating functions 



ft{x) - | 

\x, 

ott 

[x^t ■ 

X « 1, 

, X » 1 , 

(44a) 

fis(x) - | 

\x 2 , 

IX 

X « 1, 

X » 1, 

(44b) 

fl c( X ) - | 

fl, 

[o C 

X « 1, 

X » 1. 

(44c) 


In these equations, the positive phenomenological coeffi¬ 
cients a tjSjc and r] t s can be expected to be of order one. 
Note, that JF-invariance 63 dictates the prefactor (1 + 7s ) 
in Eq. (43b). 

Under the additional assumption of positivity of / 7c (x) 
it follows that (Lt x < Lt < Lt c ) 

Tc T c 

mm(a c ,l) T _ C T < | 7c (i T )| < max(a c , 1) . (45) 

Hence, the relative difference between T x and T c is small 
as 


Tx~T c 

T c 


Gi x = 


7C(3) 

7r 3 ff.Y 


« 1 . 


(46) 


Here we use the symbol Gix for the Ginzburg-Levanyuk 
number, Eq. (9), evaluated with the normal state con¬ 
ductance g x - g(L Tx ) = 2/[iTt(LT x )] » 1 . 


T c -T 
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We can further use that ft(x) and f la (x) have a max¬ 
imum ~ 1. Therefore we can estimate 


dint / t max 

d( 7 " 1 ) ~~ 2 min(a c , 1 ) ’ 

which yields 


t( L T c ) ~ t(L Tx ) 
t{Lr 0 ) 


< Gi x . 


(47) 


(48) 


In a similar way find 

7 1 ry /max 

CL In ^(jj < J'y s 

dire 1 ) - t(L Tx ) min(a c , 1 ) ’ 

and therefore 


(49) 


Z u ( Lt c ) - Z u ( Lt x ) 

z.(L T J 


<G7 


v- 


(50) 


All in all, we have determined the critical temperature 
of the fermionic system T c ~ T x with accuracy of Gi x to 
be the following function of bare parameters 


T x 

Tbcs 


= exp -2 y x 



(51) 


Here y x is the running RG logarithmic scale y at which 
the system reaches the condition (41). 


D. Ginzburg-Landau functional 

In Sec. IIC we discussed the range of length scales 
L s [Z,Lt c ] and found that the normal-state NLcrM RG 
can be applied to determine the critical temperature T c 
up to an uncertainty ~ Gi x « 1 associated with the 
critical Ginzburg region. Now we are going to ascend to 
the next level in our hierarchy of length scales and derive 
the GL functional. 

While Secs. IID 1 and IID 2 are devoted to the deriva¬ 
tion of the standard GL functional (to all orders in the 
order parameter field), in Sec. IID3 we go beyond the 
standard paradigm and derive mesoscopic fluctuations 69 
of the parameters in the GL expansion. The perturba¬ 
tive renormalization (Sec. IID 4) of the GL functional 
incorporates the effect of thermal order-parameter fluc¬ 
tuations in our theory. 


field A a (x). After this, we evaluate the resulting NLcrM 
action at q = A, which yields the sought GL functional. 

The frequency term and the quadratic Hubbard- 
Stratonovich term generate the exact GL potential 


s: 


GL 

pot 


X 


-2Z w tr[ge] 


47, 


EIa q 


J q=A 


= -27, 


+ 


m=o ■ ^7 cT 

2 hi (T/T c ) 


|A 0 


7T7 c T 

E/{s Ex/iTN 

4 

, - (-l) i (2Q!(2 2; - 1 -l) . /o , J 

§ (2Z - l)(Z!) 2 4 z 2 2i " 1 ^ ^ “ J ) 


We introduced the reduced order parameter field 6 a = 
A q /( 7 tT) and used ingredient (ii). At the equality sign 
with a dot = the divergent Matsubara sum was cut at 
e n = T x . In a similar way, we evaluate the kinetic term 
at the saddle point level: 

Sut = ^Tr[(V[[/ A (a ; )t (? [/ A ( a; )]) 2 ] g=A 

- 

[tZ?:|A ct | 2 ] 2 [-|A Q ,| + rsinh(|A a |/r)] ) 
8|A a |3 T 2[i + C osh(|A a |/T)] J A ' J j 


In view of the fact, that we are interested in tempera¬ 
tures parametrically close to the mean field temperature 
T ~ Tbkt ~ Tmfi we will keep only leading terms in the 
GL expansion: 


T 

5 = — = 
T 

^Zf x A |A Q | 2 + ||A Q | 4 + G|VA Q | 2 , 

(53) 

where 

4 Zu T 

(54a) 


A = —In—, 

7T T c 


B - 7 7C(3) 
w 2 t r 3 T 2 ’ 

(54b) 


C = 

32 T 

(54c) 

If renormalizations are discarded, we have g 

= 9d 


and Z u = 7 w/ 4 , and the parameters of the GL func¬ 
tional take the standard form , 9,59 A = t/ln(T/T c ), B = 
7C(3W(8tt 2 T 2 ) and C = ttvD(1)/(8T). 


1. The intrareplica GL functional 


2. Fluctuation corrections: Intrareplica terms 


This derivation of the GL functional for the static 
order-parameter field A Q (a:) is based on the normal state 
NLcrM at scale Lt x using (i) the local solution of the 
gap equation, Eq. (21), and (ii) the implicit definition of 
the critical temperature T c via 7 X(fr x ) = In (T c /T x ). 70 
In view of (i), we perform the same rotation, Eq. (24), 
as above, but with spatially dependend order parameter 


In Sec. IID 1 the GL functional was derived using the 
local saddle point solution of the NLcrM. However, more 
precisely, the GL functional is determined by 


S ,gl [A q ] := - In 



VQ n 2?[A Qj „,A a ,„]e- SHs[Q ’ A] , 

ntO,cx. 

( 55 ) 
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where Shs[Q>A] was introduced in Eq. (18) and 
Aa,n=o( x ) = A a(x). The functional integral is domi¬ 
nated by the local saddle point configurations leading to 
Eq. (52), while corrections are generated by fluctuations 
around these mean-field solutions. Generally, such terms 
lead to a renormalization of GL parameters A, B, C in 
Eq. (54). These effects are already taken into account in 
our approach: the solution of exact RG equations (38) in¬ 
cludes all corrections stemming from scales L e [1,Lt c ]. 
As far as corrections from larger scales are concerned, 
these will be discussed in Sec. IID 4 in the framework of 
the perturbative treatment of the GL functional. We will 
discard, however, weak (anti-) localization corrections to 
C (stemming from scales between Lt c » Lt and the de¬ 
phasing length L^). Indeed, these corrections are of the 
order 


rpWl 

—-q - Gi x In Gi x (56) 

and turn out to be subleading with respect to corrections 
discussed in Sec. IID 4. 


3. Fluctuation corrections: Disorder terms 


As shown in Appendix A, the fluctuation determinant 
associated to linear deviations from the local saddle point 
solution further leads to qualitatively new terms in the 
GL functional: these are interreplica interaction terms. 
As we explain below, most of them can be interpreted 
as random fluctuations of GL parameters A and C. The 
leading terms in the expansion |A|/T and under the as¬ 
sumption that fluctuations of A are smooth on the scale 
Lt are the following: 




■Mis 


l 

2T 2 


E f {((A{x)A{x')m a \l\^\l 

ap Jx ’ x ( 


+ 2((A{x)C{x')))\dAc\l\^\l 

+(diA* a [s x ]ii'd i 'A a ) x (d j Ap[s x ] jr d j 'Ap) X ' 

+(diA* a [s z ]ii'd i -A a ) x (d j Ap[s z ] jj ’d j 'Ap) x :l^ 

(57) 


The matrices s XtZ are Pauli matrices in the space of spa¬ 
tial coordinates i,i',j,j r e {x,y}. The correlations func¬ 
tions of GL parameters entering Eq. (57) are given by 


A(x) 

2 ~jfC(x) 


= £ ( 
D(2tt) 3 T \ 



1AA -^AC\ 6(X _ X , ) 
~1AC ICC ) 


( 58 ) 


with 



7C(3) 

(59a) 

7 AA = 

47T ’ 

7 AC = 

64’ 

(59b) 

ICC = 

7^ 2 C(3)-62C(5) 

8tt 

(59c) 


We recall that the diffusion coefficient (at scale Lt x ) is 
D = g/WZ^. Further, the parameter /3 takes values /3 = 4 
in the presence of spin-rotation invariance and /3 = 1 in 
the absence of the latter. 

Here we keep only the Gaussian white noise part of 
the distribution functions for random quantities A and 
C. The omission of interaction corrections 69 beyond the 
renormalization of D and of higher moments is justified 
in the limit Gix « 1. A more detailed study is relegated 
to a future publication. In the given approximation, our 
result for random fluctuation of GL parameter A is con¬ 
sistent with the result presented in Ref. 69. Note that 
fluctuations of A can be both interpreted as fluctuations 
of T c and of the density of states. The spatial delta func¬ 
tion is smoothened on the scale Lt, and in Appendix A 
we present the exact correlation function in momentum 
space. 

The GL functional with disorder-induced corrections 
also contains random fluctuations of higher angular har¬ 
monics, see the last two lines of Eq. (57). Since we as¬ 
sumed spin singlet superconductivity from the outset, the 
leading corrections are of d -wave type. Their distribution 
function has zero mean while the Gaussian correlation 
happens to be half as large as the one for the GL param¬ 
eter C. 


4- Perturbative renormalization 

We have now completed the derivation of the disor¬ 
dered GL functional, which is the appropriate theory at 
length scales larger than Lt c ■ We turn now to the role 
of fluctuations on scales L € [ Lt,, , £]. To investigate their 
effect, we will employ a self-consistent perturbative treat¬ 
ment. 

We observe that the mean free path of free Cooper 
pairs scattering off fluctuations of the transition temper¬ 
ature is 


^bosons ~ £\tc\g » ((, 


where the reduced temperature is 


e c 


T-T c 

T c 


(60) 


(61) 


Therefore, in the considered regime of length scales, the 
Cooper-pair kinetics is ballistic. 

Before turning to self-consistency, we present the sim¬ 
ple perturbative correction to the mass term 
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A 


A 


1 + 2 —A In 



with 


A + 2 A— In 

e c 



A = 


4 


T-T c 

T r 


+ 2 Gix In 


(£)] 


(63) 


This correction stems from quartic term of the clean GL 
functional, Eq. (53). The analogous diagram from the 
quartic interreplica term, Eq. (57), is smaller by addi¬ 
tional factor of Gix■ We note that Gix/e-c is the same 
as the parameter 'y c (LT x )t(Lx x ) from the analysis of 
NLcrM RG, Sec. IIC. In the second line we disentangled 
effects of the two dimensionless parameters in the theory 
(Gix/tc and Gix)- The effect of fluctuations is a reduc¬ 
tion of the transition temperature, which turns out to 
be 9 


e » Gix- As explained in Ref. 3, in this case the RG has 
to be stopped at the scale Lt- The corrections to the 
kinetic GL parameter C -*• C can be associated to the 
conductance 


9 


9 



( 68 ) 


Here, Gip is the Ginzburg-Levanyuk number with g(Lr) 
evaluated at scale L T . The expression (68) can be inter¬ 
preted as the (physical) conductance above Tmf and is 
of the same form as was presented in Eq. (45) of Ref. 3, 
provided Ii^L^/Lt) —ln(Gix/|e|). 

Furthermore, we conclude that the GL coherence 
length, Eq. (11c), is more accurately given by 


£gl (T) 


m D(L T ) 
\| D(L Tr y 


e » Git- 


(69) 


T MF = T c (l-2Gi x \lnGi x \). (64) 

Further, the shift of transition temperature is taken into 
account self-consistently. This results in replacing e c -*■ e, 
where 


e=(T-T M F)/T M F, (65) 


and A -*■ A in the last term of Eq. (62). Summing up all 
of these effects leads to a replacement 


A 



( 66 ) 


in the GL functional, Eq. (53). The factor in square 
brackets in Eq. (66) can be absorbed into a redefinition 
of the order parameter field A a . Then, the self-consistent 
perturbative treatment yields a GL functional of the form 
of Eq. (53) and Eq. (57) with renormalized coefficients 


The zero temperature coherence length is £(0) ~ 

yD(L T JlT c . 


E. Disordered 0(2) model, superconducting 
density and vortex unbinding transition 


In the previous section we derived the effective GL 
functional and incorporated the effect of fluctuations 
from scales L e [ Lr c , £]. We now reach the largest scale 
at which a crossover of theories takes place, namely the 


GL coherence length £ = yjC/(2\A\). In order to obtain 
an effective description in terms of the U(l) NLcrM simi¬ 
lar to Eq. (3), we need to determine the mean-field order 
parameter. 


1. Order parameter fluctuations 


A -+ 
B -» 

C -+ 


I=^e, 

7r 

B = - 


B 


C = 


1 + 2^ In 

C 


m 


1 + 2 In 


(ffir 


(67a) 

(67b) 

(67c) 


The parameteres JaAiIac an< 4 lee are renormalized 
analogously to the coefficient B. 


5. Metallic side of the mean field transition 

We note that the procedure employed up to now is 
equally applicable for T > Tmf and T < Tmf- Let us 
briefly concentrate on the region of weak fluctuations on 
the metallic side of the superconducting transition, i.e. 


The mean-field equation for the order parameter field 
is 

-V(CVA M f(®)) + AAmf(®) + S|A M F(®)| 2 A M F(a;) = 0. 

(70) 

Here, the quantities C(x) = C_ + SG(x) and A(x) = 
A + 6A(x) fluctuate randomly according to Eqs. (58) 
and (67). Iterative, perturbative solution of this equa¬ 
tion leads to Amf = Amf + Wmf(*) 


|A 


MF - 


N 


= 27tX , mf -< 
±3 


, 2 |e| 

M 7C(3) . 


1 + 2 —A i n 


and 


(71a) 


^A MF (a:) = -A mf 



6A(x')e ik ( x - x ' ) 
C[k 2 + yi] 


(71b) 
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Therefore, the order parameter field fluctuates 69, ' 1 on the 
scale of the coherence length 


({ Amf (3:) Amf (*£ ^))) ((^MF(®)A^[ F (a:'))) 


A 2 

U MF 


|A 


MF 


((^mf( 3: )^mf( 3; '))) 

[AmfI 2 


= ICa(x - x') 


(72) 


with Gaussian correlation function 


ICa(x) = _ 2 


PlAA 


L 


Akx 


C D(2tt) 3 T Jk [k 2 + £~ 2 ] 


PlAA C 2 


C D(2tt) 3 T 


47r 


,(M). 


£ 


f3 jaa tt 6 
49C 2 (3) L 


Gi x 


i 2 


D5(x - x') 
27tT 


(73) 


The fluctuations of K are given by 

SK(x) _ SC(x) + iIAmf)^) + <5A^ F (a;) 
K (7 A M f A^ f 


(76) 


The fluctuations of K are goverened by order parame¬ 
ter fluctuations: 


(6K(x)SK( x')) 


PlAA 

4 


7 r 3 2 D6(x - x') 
7C(3) J 2 ttT 


(77) 


Note that the delta function is smoothened on the scale 
£. Mesoscopic fluctuations of the GL parameter C lead 
to subdominant contributions to (5K(x)6K(x')). We 
explicitly checked, that terms proportional to 7 ac or "/cc 
are negligble as compared to those proportional to 7 aa 
in view of additional small factors of |e| and/or Gix- 


where K t is the first modified Bessel function. The white- 
noise approximation in the last line of Eq. (73) is justified 
when the physics on spatial scales much larger than £ is 
considered. 

It should be noted, that the inclusion of higher Matsub- 
ara harmonics can lead to long range ln 2 (fc£) tails in the 
Fourier transform of K.a(x).' 1 Such effects are beyond 
our accuracy in the determination of Tbkt+ see the dis¬ 
cussion in Sec. IIF. Further, we explicitly checked that 
fluctuations of C{x) do not contribute to fluctuations of 
Amf(*)- 


2. Mean field stiffness 


3. BKT transition temperature 

We remind the reader that the clean U(l) NLctM sup¬ 
plemented with topological excitations (vortices) exhibits 
a phase transition which is driven by the logarithmic in¬ 
tervortex interaction. The latter is encoded in the follow¬ 
ing RG equations for the mean stiffness K and fugacity 

%V 1 


dK- 1 

dy 

dz v 

dy 


2(1 -K)z v , 


(78a) 

(78b) 


Finally, we find the mean field stiffness, which fluctu¬ 
ates on the scale of £, to be 


K(x) C(x) 2 K + SK(x) 

- = rr, IAmFWI “ -• 

IT 1 IT 


(74) 


Following the notation from above, K denotes the aver¬ 
age stiffness and SK(x) the small fluctuations. Replacing 
C and Amf by their respective expectation values leads 
to 


K 

7T 


3 I Amf I 
32 T 2 

47t Gix 


1 + 2 


Gix 


1 + 2—A in 



(75) 


Here we have introduced |A(J ( f | 2 = -A/B 

-e8n 2 T' 2 /7£(3), the mean field gap without corrections 
stemming from length scales L e [i Tc , £]. 

Thus, the mean superconducting density K(e) (e < 0) 
is determined by the normal state conductivity at the 
same distance from the transition (i.e. g(-e) obtained 
by reflecting e = (T - Tmf)/7~mf about the origin), see 
Eq. ( 68 ). 


where y = In (L/t;). The fixed point defining the vortex 
unbinding transition is given by ( K* ,z*) = ( 1 , 0 ). 

Generally, the effect of disorder on the vortex physics 
in the U(l) NLtrM is twofold: first, the random stiff¬ 
ness modifies the interaction between vortices and rota¬ 
tionless supercurrents as well as between vortices them¬ 
selves. Second, since the stiffness determines the core 
energy S core = - In z v = a z K , fluctuations of the stiffness 
can lead to pinning of vortices. The coefficient a z is a 
number of order unity in realistic superconductors, 3 ' and 
the fugacity z v acquires a log-normal distribution. 

In the present case, the effect of disorder on the BKT 
transition temperature is relatively weak: the quartic in¬ 
terreplica terms are RG irrelevant. In addition, since 
K(x) fluctuates on the scale £, vortex pinning is negli¬ 
gible for the following reason: The size of a vortex core 
is of the same scale £ and by consequence fluctuations 
largely cancel out. 

Technically, the effect of disorder can be taken into ac¬ 
count by partial contraction of the quartic interreplica 
terms . 72 It turns out that intrareplica quartic terms 
(thermal fluctuations from length scales L > £) dominate 
upon disorder effects by a factor of Gi 3 /. 

Using this observation and simple rescaling 
arguments 73 in the formal definition of the bare 
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stiffness we find that Rb are depends only on the single is conserved under the BKT RG, Eqs. (78). On the nor- 
parameter Gix/ |e| mal conducting side, the limit tbkt -*■ 0 implies 


/W 1 |e| iGix 

7 T 47 t Gix K \ |e| 


(79a) 


with fx(x)/x being a monotonically decreasing function 
and 


fi< (a;) 


l-2xln:r, fora;«l, 

< x VK ( t\k < 1), for x » 1 . 


(79b) 


In the limit z v -*■ 0, the critical stiffness is K * -»■ 1. The 
transition temperature Tbkt is governed by the solution 
«bkt of /a'(1/o!bkt) = 4/q!bkt, he. 


|cbkt| 


Tmf ~ Tbkt 
Tmf 


a BKT G'ix- 


(80a) 


Even though the precise value of q;bkt remains unknown, 
the asymptotic terms in the limit |e| » Gix are sufficient 
to identify «bkt as a number of order unity. 

In a realistic superconductor close to the transition, the 
fugacity z v ~ 0(1) is finite and thus the system under¬ 
goes the phase transition far from the critical end point 
(K*,z*) = (1,0). This does not alter, however, the result 
(80a) with «bkt understood as an unknown coefficient 
of order unity. For completeness we repeat that 


Gix = 

7<<3) «i 

(80b) 



Tmf = 

T c (l-2Gi x \\nGi x \), 

(80c) 


with gx being the NLctM coupling constant evaluated 
close to the superconducting instability of the diffusive 
Fermi liquid occurring at the energy scale T c . The ex¬ 
pression for the vortex unbinding temperature, defined 
via Eq. (80a), is one of the major results of this paper. 

We speculate now on the behavior of the vortex un¬ 
binding temperature, when the system is tuned close to 
the quantum critical point [characterized by ( 7 c (Tt’ c = 
oo ) 1 t{LT c = °°)) = (—oo, tgiT ~ 1)] of the NLctM RG from 
Sec. IIC 4. It appears very plausible that Eq. (80a) de¬ 
scribes the first term of a geometric series, such that 


Tbkt 


Tmf 

1 + a B KTGkY 


(81) 


A result of this form was first obtained by Beasley, 
Mooij and Orlando, 26 who however omitted all fluctu¬ 
ation and disorder renormalizations. (In their approxi¬ 
mation Tmf = TbcSj <tbkt = 4, and Gix - Gi.) 


4- Correlation length 

It is easy to check that the quantity 

c(I<) = 4 [K- 1 - 1 - In AT -1 ] - z 2 v (82) 


/ tbkt 
V 26/ 2 Gix 


(83) 


We remind the reader of the definition tbkt = (T - 
Tbkt)/Tbkt- The fitting parameter b is chosen such that 
the solution of RG Eqs. (78) with c(A') given by Eq. (83) 
leads to the same form as Eq. (lib) 


£bkt(T) ~ £(T B KT)e ^ TB ' KT , t B kt « Gix (84) 

for the correlation length.' 4 As was discussed in detail 
in Ref. 37, the parameter Gi (or better Gix) naturally 
appears in the standard definition of both BKT and GL 
coherence lengths. 


5. Temperature-dependent stiffness 

We now focus on the superconducting side of the tran¬ 
sition and discuss the temperature dependent mean stiff¬ 
ness. We consider three different temperature regimes: 
(i) T « Tbkt, (ii) Gix « |tbkt| « 1 and (iii) |tbkt| « 
Gix- 

In the immediate vicinity of the transition [regime 
(iii)], Eq. (83) implies 


A"(tbkt) 


1 

1 - fb\/\ T BKT\/Gix 


(85) 


In contrast, in regime (ii), renormalization effects are neg¬ 
ligible and Eq. (79) describes the temperature dependent 
stiffness. 

Finally, we consider the regime (i) of lowest tempera¬ 
tures, T « Tbkt- It is useful to first clarify the hierarchy 
of length scales occurring in this regime. As compared to 
the case -tbkt « 1 depicted in Fig. 1, now the following 
scales are almost identical: Lt c ~ Ta ~ £■ Therefore, the 
large window Lp c « L « £ of dominant thermal fluctu¬ 
ations disappears and the stiffness is determined by the 
outcome of the NLctM RG, 


AT 

7T 


felA tanh A. 

8 2 T 2 T 


( 86 ) 


Within our theory, A = A(T) » T is the temperature- 
dependent spectral gap, which has standard BCS form, 
but with renormalized T c , as dictated by the NLctM RG. 
This expression incorporates the quantum corrections 
stemming from scales L < Lp c - The low-temperature 
stiffness (superconducting density) including dominant 
quantum corrections to T c was recently determined in 
Ref. 75 using an alternative approach (self-consistent di¬ 
agrammatic technique). 
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F. Summary and discussion 


A. Interpolating function for the resistance 


The final result for the vortex unbinding temperature 
is visualized in Fig. 2. For concreteness, we consider a 
system with long-range Coulomb interaction and strong 
spin-orbit coupling. The transition temperature is dis¬ 
played as a function of the Drude resistance tp. Devia¬ 
tion from the horizontal line, T/Tbcs = 1, implies a viola¬ 
tion of Anderson’s theorem. The three curves correspond 
to Finkelstein’s solution, Eq. (6) (blue dashed line), to 
the temperature Tx of the superconducting instability, 
Eq. (51), as resulting from the full one-loop fermionic 
RG (red clot-dashed line), and to the vortex unbinding 
temperature Tbkt, Eq. (81) (green full line). We see 
that there is a substantial difference between the three 
curves. The relative difference between Tx and Tbkt 
is of the order of Gix\ lnGt A |, see Eqs. (46), (64), and 
(80a). This difference can serve as a rough estimate for 
the width of the BKT fluctuation regime in the temper¬ 
ature dependence of resistivity, which is of order Gix , as 
will be discussed in Sec. III. 

Let us emphasise once more that, as seen in Fig. 2, the 
found mean-field transition temperature Tmf, and conse¬ 
quently the true transition temperature Tbkt, can differ 
strongly from the BCS temperature Tbcs- The domi¬ 
nant effects behind this difference are those incorporated 
in the NLctM RG. In addition to them, there is a contri¬ 
bution of GL fluctuation corrections. The randomness of 
the stiffness inherited from “mesoscopic” fluctuations of 
GL parameters, Eq. (59), is parametrically smaller and 
thus of minor importance for the finite-temperature su¬ 
perconducting transition that we are exploring. This ap¬ 
plies also to ln 2 (fc£) tails' 1 in the correlation function of 
K(x), Sec. HE 1, which are subleading and do not affect 
Eq. (80a). 

Therefore, our theory implies that in homogeneously 
disordered films the long-wavelength theory describing 
the finite temperature vortex unbinding transition is a 
theory of clean bosons. In this paper we showed how 
this theory emerges from the underlying theory of dis¬ 
ordered fermions. Our derivation of parameters of the 
effective U(l) theory takes into account all essential dis¬ 
order, quantum-interference, and interaction effects. 

It is worth pointing out that we have assumed a short- 
range disorder. Spatial variations of the chemical poten¬ 
tial and/or the BCS-coupling constant with a character¬ 
istic length large compared to £(Tbkt) can be incorpo¬ 
rated in the resulting bosonic theory, which will yield a 
disordered U(l) NLcrM. This will lead to an additional 
suppression of Tbkt- 3 ' 


The asymptotic behavior of the resistance in the vicin¬ 
ity of the BKT transition is given by generalization of 
the result obtained in Ref. 36. 

P ( T ) - -KH x | ) ’ T BKT«Gi A - (87) 


and by Eq. (68) 


P(T) 


t(L T ) 


1 + 4 


£hn(T) \ 2 ln / Chn(T) \ 
£(Tbkt)/ \£(Tbkt)/ 


tbkt » Git- (88) 


Here £hn is a generalization of the Halperin-Nelson 
length, 36 cf. Eq. (11a), 


{„»(!■)- ffifsWi,AG, (89) 

b \ tbkt ) 

inasmuch it contains all renormlization effects discussed 
in this work. We introduced an interpolating function for 
the Ginzburg number Gir = 7()(3)t(T)/(27r 2 ) by means 
of t{T ) given by f(T) = f(T-r) for T > Tx and by the 
extrapolation 

t(T) = Ta-Gi a 9t^(Tt)|t x tanh 

for temperatures below Tx- 

The vortex-generated resistance in the immediate 
vicinity of Tbkt, Eq. (87), contains an unknown pref¬ 
actor A” 1 of order unity. The resistance far from the 
superconducting transition includes the fluctuation cor¬ 
rections, Eq. (68), which are essentially the anomalous 
Maki-Thompson correction. In 2D, they dominate upon 
Aslamazov-Larkin corrections. We propose the following 
interpolating function for the resistance which reproduces 
correctly both limits (87) and (88): 


T-T X 
TxGix J 


+ t(L Tx ) (90) 


P(T) 


t(T ) 


1 + 2 


( ;hn(t) y 

V s(Tbkt) / 


In 


„ A „/2 



(91) 


Let us comment on the width of the crossover region 
from metal to superconductor in the temperature depen¬ 
dence of resistance. As is indicated by Eqs. (87) and (88), 
our formalism predicts that this crossover happens within 
a window of relative size Gix around the mean field tem¬ 
perature Tmf- The only temperature scale arising in our 
analysis of the crossover is |e| ~ Gix- 


III. TEMPERATURE DEPENDENCE OF 
RESISTANCE 


B. Electrical vs. heat transport 


This section is devoted to the temperature dependence 
of the resistance in a 2D metallic film at the verge of 
superconductivity. 


The final result for the resistivity, Eq. (91), is displayed 
in Fig. 3. Again we focus on the case of a system with 
strong spin-orbit coupling and Coulomb interaction. 
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T / T BCS T / T BCS 




FIG. 2: Transition temperature for a homogeneously disordered superconducting film with Coulomb interaction and strong 
spin-orbit scattering as a function of the Drude resistance to = 1/(vr pr). On the left (right), the mean free time r is varied 
(fixed) while the chemical potential is kept fixed (is varied). Blue dashed: Finkelstein’s approximate solution for T c , Eq. (6). 
Red dotted: Temperature Tx, Eq. (51), of the instability of the fermionic system. Green solid: Vortex unbinding temperature, 
Eq. (81). Green and red curves are obtained by numerical solution of RG equations. The only free parameter in the left 
(right) plot is Tbcs/m = 0-01 [In(TBcsT) = -0.2]. The origin of the different behavior in left and right plots is due to the strong 
dependence of 7c(0 on to in the left graph. 


When the transition is approached, the physical resis¬ 
tance p starts to deviate from the NLerM running cou¬ 
pling constant t and eventually rapidly decreases to zero. 
The difference between the two quantities is due to inelas¬ 
tic processes, in particular fluctuation corrections, which 
are beyond the RG scheme of the NLerM. The quantity 
t(L ) serves as a starting point for the calculation of var¬ 
ious physical observables subjected to these corrections, 
including the electrical resistance, Eq. (91), and the su¬ 
perconducting density, Eq. (79). We interpret the NLerM 
coupling constant t(L) [or t(T)] as the electronic contri¬ 
bution to the thermal resistance of the system. Indeed, it 
is known that the quasiparticle contribution to the ther¬ 
mal transport coefficient n has no singular corrections 
at the onset of the transition. 79 Technically, this is re¬ 
lated to the extra frequencies at heat-current vertices 
(see Appendix B), each producing an extra factor y" 1 
in the inelastic contribution to the thermal conductivity 
as compared to the electrical conductivity (see Ref. 3 for 
the NLerM calculation of the latter). In particular, the 
thermal conductivity k can be readily derived within our 
NLerM formalism at the saddle-point level (see Appendix 
B for the technical details), with the result 


g_T r°° t 2 d( 

4t r Ja/t cosh 2 (C/2) ’ 


(92) 


and is thus related to renormalized g, and consequently, 
to renormalized t. 

To summarize, well above the transition both electric 
resistivity and thermal resistivity t follow the NLerM RG. 
When the transition is approached, the electric resistivity 
p gets strongly suppressed due to fluctuation effects and 
eventually becomes zero at Xbkt- On the other hand, 
the NLerM coupling t, and thus the thermal resistivity, 
remains finite in the vicinity of the transition tempera¬ 
ture Xbkt- 


tip 



FIG. 3: Temperature dependent resistance of a homoge¬ 
neously disordered metallic film close to the superconducting 
transition. Here, the case of strong spin-orbit coupling and 
Coulomb interaction is considered. Blue dashed curve: tem¬ 
perature dependence of the NLerM coupling constant f, which 
determines the prefactor of Eq. (92) for the thermal conduc¬ 
tivity. Green solid curve: temperature dependence of the 
physical resistance. The BCS transition temperature Xbcs is 
marked by the black dashed vertical line, while the mean field 
transition temperature Tmf is marked by a dot-dashed line 
surrounded by a gray Ginzburg-Levanyuk window. Note the 
logarithmic scale of the temperature axis. This plot was gen¬ 
erated using Eq. (91) with the bare values ln(TBcsT) = -0.2 
and tn = 0.1, as well as phenomenological constants obkt = 4, 
A v = 7r 4 /56((3) and 6=1. 


It should be mentioned that in the model of diffusive 
electronic system with long-range Coulomb interactions 
there are additional logarithmic corrections to the ther¬ 
mal conductance k. 80 These corrections, however, are 
nonuniversal and beyond the RG scheme, 81 so that they 
are not accounted for explicitly in our treatment. 
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C. Finite size effects and the role of 
electromagnetic fields 

Finite size effects may lead to a smearing of the BKT 
transition. First, as for any phase transition, the lin¬ 
ear dimension L a of the film leads to a natural infrared 
cut-off. Second, a finite size effect which is specific to su¬ 
perconducting films is due to the finite thickness d + 0, as 
has been already mentioned in Sec. (IIC 1). Due to the 
interplay with electromagnetic fields, a finite thickness 
leads to an infrared cut-off for the logarithmic vortex- 
antivortex interaction 38 at the length scale of A 2 /d, where 
A is the London penetration depth. In the context of the 
present treatment, this means that RG has to be stopped 
at the smaller of the two scales {A 2 /d,L n }, which at the 
same time is the maximal value for £bkt- The conse¬ 
quence is a rounding 37 of the sharp resistance drop near 
T B kt in Fig. 3. 


IV. SUMMARY AND OUTLOOK 

To summarize, we have developed a theory for the 
vortex unbinding transition in homogeneously disordered 
superconducting films. This theory incorporates the ef¬ 
fects of quantum, mesoscopic and thermal fluctuations 
stemming from length scales ranging from the super¬ 
conducting coherence length down to the Fermi wave¬ 
length. Using the developed theory, we determine the 
dependence of essential observables (including the vortex 
unbinding temperature, the superconducting density, as 
well as the temperature-dependent resistivity and ther¬ 
mal conductivity) on microscopic characteristics such as 
the disorder-induced scattering rate and bare interaction 
couplings. More specifically, our key results are as fol¬ 
lows: 

1. We have performed a consecutive mapping 
of emerging effective field theories, starting 
from the Fermi-liquid theory at a short scale 
through the normal-state interacting sigma model, 
superconducting-state sigma model, and Ginzburg- 
Landau theory, to the U(l) theory of supercon¬ 
ducting phase fluctuations. 

2. As a result of this procedure, we obtained the su¬ 
perconducting density in the temperature regimes 
(i) T « Tbkt, Eq. (86), and (ii) Gix « |tbkt| = 
\T - TbktI/Tbkt « 1, Eqs. (75), (77), (79), which 
include effects from quantum, mesoscopic and ther¬ 
mal fluctuations. We find that in regime (ii) the 
superconducting density at temperature T is de¬ 
termined by the the normal state conductance at 
the temperature 2Tmf - T obtained by reflecting 
about the mean field transition temperature Tmf- 
The scale Lt x and the associated Ginzburg num¬ 
ber Gix are provided by the fermionic nonlinear 
sigma model RG. 


3. Combining our result for the superconducting den¬ 
sity in regime (ii) with the RG treatment of the 
BKT transition, we have found the vortex unbind¬ 
ing temperature Tbkt, Eq. (81). Figure 2 depicts 
the transition temperature as a function of exper¬ 
imentally controllable parameters in the case of a 
system with strong spin-orbit coupling. Further, 
we extracted the behavior of the stiffness in the 
immediate vicinity of the transition |tbkt| « Gix , 
Eq. (85). We thus provided results for the super¬ 
conducting density in all regimes below Tbkt- 

4. We have proposed a function, Eq. (91), which inter¬ 
polates between vortex-dominated (exponentially 
small) resistance close to the vortex unbinding tem¬ 
perature and the Maki-Thompson-likc fluctuation 
resistivity representing the dominant fluctuation 
correction further from the transition. We find, 
that the resistance drop occurs in a small temper¬ 
ature window of strong fluctuations of relative size 
Gix • Our result for the resistance is plotted in 
Fig. 3 for the case of a system with strong spin- 
orbit scattering. 

5. We have identified the running charge t(L ) of the 
NLctM RG scheme of our theory with the electronic 
contribution to the thermal resistance of the sys¬ 
tem, see Eq. (92). 

6 . We have performed a derivation of mesoscopic fluc¬ 
tuations of Ginzburg-Landau coefficients, including 
both the mass term and the prefactor of the kinetic 
term, for both cases of the preserved and broken 
spin-rotation invariance. 

Before closing the paper, we briefly discuss some of 
perspective directions for future work. 

1. Our results, and in particular Eqs. (81) and (91), 
should be useful for the analysis of experimental 
data on Tbkt and on resisitivity near the transi¬ 
tion, including recent measurements 47-53 and ex¬ 
pected future experiments. It would be also very 
interesting to prove experimentally our observation 
that the distance (Tbkt -Tmf)/Tmf should be of 
the order of the dimensionless thermal resistivity. 

2. This paper was devoted to the thermodynamic 
phase transition into the superconducting state. 
While renormalizations on intermediate scales were 
of quantum origin, the physics in the close vicin¬ 
ity of the transition was governed by thermal fluc¬ 
tuations. The corresponding transition line ends 
at T = 0 at the point of superconductor-insulator 
transition, 1 which represents a prominent exam¬ 
ple of a quantum phase transition and is driven 
by quantum fluctuations. Contrary to the finite- 
temperature transition controlled by a Gaussian 
fixed point, the fixed point governing the quantum 
superconductor-insulator transition is at strong 
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coupling, which makes its controllable analytical 
exploration an extremely difficult task. While con¬ 
siderable progress has been achieved within theories 
of disordered bosons 76- ' 8 and within the fermionic 
NLctM RG, see Ref. 3, development of a control¬ 
lable theory of the critical behaviour at this quan¬ 
tum phase transition remains an outstanding chal¬ 
lenge for future research. 
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Appendix A: Derivation of disordered GL theory 
from the NLctM 


In this appendix, the interreplica corrections to the 
GL functional are derived. To this end, the functional 
integral Eq. (55) is evaluated beyond the saddle point 
approximation. 

Starting point for the derivation of the interreplica 
terms of the GL functinal is the action (18), which 
yields, after integration of Hubbard-Stratonoviclr fields 
A an with n + 0, 


Refs. 3,23). This leads to 


5 = J^Tr{(VlF) 2 + [G A Vt/i,Af 


+2 [W,VW]U A VUl 

+2VlE[t/ A VC/^,A] + 2W[[U a VUI : A],UaVU{] 
- [t/ A Vt/l, A] [UaVU^AW 2 ] + [U A VUl, wf } 


-2Z,Tr 




4Z "EiA„| A s-y- 


■KjcT ‘ 


(A2) 


We first present the strategy for the calculation of {{AA)) 
in the noninteracting case and consider interaction effects 
and other coefficients ((AC)) and (( CC)) afterwards. 


1. Fluctuations in GL coefficient A: noninteracting 

case 


To obtain the coefficient ((AA)), we can omit all VC/ A 
terms in Eq. (A2) as well as the interaction terms. Then 
the action can be diagonalized with the parametrization 


IT = 



where we expand (er,- = (l a ,d)j) 


^ / 7a/3 _ - a/3 ' 

,.,a/3 _ _ I a nin 2 ;j lC nm 2 ;j 

a ni,n 2 9 I -~a/3 ,a0 

\ to nin 2 ;i n 1 n 2 ',j , 

\ / ia(3 -aft 

/3a _ A? [ a n 1 n 2 ;j 1C nm 2 ;j 

" 2,ni 2 1 ;A a/3 d a/3 

\ nm 2 ;j u n 1 n 2 ;j 

Here, the notation A j means 



f-1, if j 6 {1,2,3}, 
11 , else. 


(A3) 


(A4a) 

(A4b) 


(A5) 


S = 


X|.r(W)= 


2Zu.tr 


Q 


e + 


4 Z„ 
7T7 C T 


S A«Aa + Sfr 




(Al) 


The collection of interaction terms, Eqs. (15c)-(15e), 
(keeping only nonstatic parts in the Cooper channel) is 
abbreviated as 5^’^”*°. 

As a next step, the space-dependent Bogoliubov ro¬ 
tation [analogous to (24)] is performed. For the deriva¬ 
tion of interreplica terms in the GL functional, the fluc¬ 
tuations around the saddle point are parametrized by 
q = IT + A\/l - IT 2 rs A + W - AW 2 /2 with IT Hermitian, 
IF A = -IT A and IT = ti 2 W T ti 2 (same notation as in 


We use the convention that negative (zero or positive) 
Matsubara indices are denoted by even (respectively, 
odd) subscripts. The diagonal action is 


*[<*.=] = tzJ E vF2(9)]I[, 

lu J a n\no. ’ 


’q mn 2 
a(3 

0<j<3 


\d a > 3 d afl 

l“ , n 1 n2-,3 Uj n 1 n 2 ;, 


where 




, ~otQ a. 

4- r C 

3 n 1 n 2 \3 


+P /3 


A m 


a/3 _2 , £ m + e n 2 


D 


(A7) 


and D - g/16^ w . For convergence reasons, the complex 
diffuson fields have the property dj = Xjd* and analo¬ 
gously for cooperons. 
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The integration of the modes yields a fluctuation deter¬ 
minant, or, equivalently, a fluctuation correction to the 
action 

6S = 2(3£ £ TVln [V^q)]^ (A8) 

a,/3 711,712 

It is worth noticing that the action (A6) is already di¬ 
agonal in the multiindex (a, /?, 711 , 77 . 2 ). Thus the trace 
in Eq. (A8) does not include trace over replicas and/or 
Matsubara indices. However, the action (A6) is not diag¬ 
onal in momentum and coordinate space and the symbol 
“TV” refers to an operator trace in this space. 

The parameter (3 is (3 = 4 in the quaternionic case with 
spin-rotation invariance, with a singlet and three triplet 
modes contributing. For broken spin-rotation invariance, 
when only the singlet mode contributes, we have (3 = 1. 

We now expand the expression + |A a (x)\ 2 

in small |A Q |/T. Recall that we are interested in in¬ 
terreplica correction to the GL functional up to quartic 
order. It is thus sufficient to keep the quadratic order 
in the expansion of (the quartic order would lead to 
quartic mtrareplica contributions with prefactor Nr -*• 0) 


£ 


a 

n 


£n| + 


|A“| 2 


(A9) 


The expansion of £“ leads to the approximate inverse 
propagator 




|A“| 2 | A' 5 1 


+ 


2 |e„i| 2 |e n 


Here, uj ni2 = e„, - e 


(A10) 

n2 and V denotes the propagator in a 


normal metal obtained from V A by setting A = 0. 

We expand the logarithm, Eq. (A8), to quartic order 
in A and keep only terms that survive the replica limit: 


SS GL 


r |A a | 2 |A^| 2 , 

£3 Jx,x' 2D 2 e ni \e n2 \ 

711,722 

*T>(x - x',u> ni2 )V(x' - X,U>n 12 ) (All) 


As will be discussed below, the evaluation of the Matsub¬ 
ara sums leads to the correlation function for fluctuations 
of the GL coefficient A. 


2. Fluctuation of GL coefficient C and interaction 
effects 

We now return to the derivation of the full, disor¬ 
dered GL functional, keeping terms up to quartic order 
in A/T. To this end, we consider Eq. (A2). It is appar¬ 
ent from the above derivation of the disorder correlations 
((A(x)A(x 1 ))) that interaction terms do not di¬ 

rectly play a role for the derivation of the disordered GL 


functional. Indeed, interaction effects only appear in the 
propagators \V~ A 1 (g)] ni ■ These terms cannot lead to 
interreplica terms, as is evident in view of Eq. (A8). Fur¬ 
ther, the first line of Eq. (A2) produces the kinetic term 
of the normal state diffusion propagator and the stan¬ 
dard (saddle point) kinetic term of the GL functional. 
The terms of the third line of Eq. (A2) are unimpor¬ 
tant, since they only affect the replica diagonal channel. 
Therefore, in addition to terms from TriAW 2 we need to 
consider the following terms, of which the first two (last 
two) stem from line 2 (respectively, line 4) of Eq. (A2): 

5S[W, A] = J^Tr{2[W,VlT][/ A V[4 

-[U A VUi,A][U A VUi,AW 2 ] 

+ [U A VUl,W ] 2 } (A12) 

As before, we expand to the second order in A/T. This 
yields 

U A VUl = -i £ IxpwVA“ + /“t 2 1 . (A13) 

n> 0,ct { Z \ e n\ ) 

Here we have introduced the following matrix in Narnbu 
space: 


A“ 


0 -iA a \ 

iA a ’* 0 ) ’ 


and the vector field 

fn ~ y-y (A“VA“’* - A“’* VA“). 

8 e„ 


(A14) 


(A15) 


We first consider the second and third term in the square 
bracket of Eq. (A12). It is sufficient to keep terms up to 
linear order 0( A*) in U A S/U'^. We thus obtain from the 
second and third terms 


^2,3 = £ Tr 


+ - Y Ti¬ 

le L -‘, 

a, a. 
n,n>0 


a;n> 0 


Pi° 


\ 2 |Cn| / 


2|e 

-iA. 


W 






(A16) 


The effect of these terms is twofold: The first term 
from £> 62,3 yields a shift of the propagator (A7) 

f P ^[V-\q)] aP 

J 721,722 L 7 J 721,722 


s [^(*)C 2 4 


«/3 r«-i,-xW> 1 /|VA“| 2 , |VA^| S 


) 


(A17) 
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Second, the last terms from < 5 <S 2,3 produce couplings be¬ 
tween cooperons and diffusons of opposite Matsubara fre¬ 
quencies: 
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a,/3 l VA*’ Q VA *’ /3 0 


), (A19a) 

j. (A19b) 


We now consider the first term (with the single com¬ 
mutator) in the square bracket of Eq. (A12). Since 
IT 2 and W\7W are block diagonal in Matsubara space, 
TrAj,W 2 = 0 = TrAyWVWA Thus only terms stemming 
from the term proportional to /“ in survive the 

trace operation. The first term thus yields the contribu¬ 
tion 


f p g {Cfn 2 ;, [/£(«) - /£,(«)] ( 2 P + 

+5 “fn 2 ;j [/“i(®) + /&(«)] (2p + 9) c “fn 2 ;j} 


(A20) 


In the replica limit we will be only interested in contri¬ 
butions from the crossed terms which cancel out between 
cooperons and diffusons in view of the different sign of 
contributions. We therefore do not need to consider SSi 
any longer. 

The fluctuation determinant leads thus to a correction 
of GL action of the form 


SS 1 
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Here, the first two lines stem from diffusons and the third 
and fourth line from cooperons. The symbol Tr denotes 
trace in momentum/coordinate space and in the space of 
the 2x2 matrices introduced in Eq. (A18). It does not 
include the trace over replica and/or Matsubara indices. 
The expansion of the trace to fourth order in A leads to 

SSdi = - E f V (x - x',u ni2 )V(x' - x,co ni2 ) 

a,/3 Jx ’ x ' 
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+ (va q va / 3 ) x (va* va;)^]J . 

(A22) 

This expression is the origin of Eq. (57) of the main text. 


3. Discussion of disordered GL functional 


All in all, our derivation yields the following inter¬ 
replica corrections to the GL functional to the leading 
order in the expansion |A|/T: 
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(A23) 


We will be interested in the limit when the field A is 
smooth on the length scale of Lt- In Eq. (57) of the 
main text we presented only the delta correlated part of 
fluctuating GL parameters. Here we will go beyond this 
approximation. In general, the Gaussian statistics of GL 
parameters is dictated by mean values given in Eqs. (54) 
and fluctuations 
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It is convenient to present the correlation functions en- 
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tering Eq. (A24) in Fourier space: 


Appendix B: Derivation of the thermal conductivity 
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In the approximate evaluation, we only kept terms yield¬ 
ing terms d n (A/T) m with n,< m < 4 in the GL func¬ 
tional. We will use the notation jaa for 7 aa(q - 0) 
and analogously for 'yac, ICC- The sums entering these 
expressions are 
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The eigenvalues of the matrix of 7 are positive in the 
limit Dq 2 /T « 1. At q = 0 they are approximately 15/47r 
and 2 . 8 / 47 r. In the main text, we further used 


(d i A* a d i A /3 )(dj A a dj 
+ (d i A a d i A f} )(d j A* a d j A* fl ) = 
(d i A* a d i A a )(d j A f3 d j A* 0 ) 

+ (di A* [s x ]ii'Si' A a )(dj Ap [s x ] jj-dj- A*p ) 

+ (d i Al[s z ]i i 'di-A a )(d j A 0 [s z ] jj 'd j 'Ap). (A27) 


In this section, we present the derivation of the quasi¬ 
particle contribution to thermal conductance in the dis¬ 
ordered superconductor. Upon taking A -> 0, this dis¬ 
cussion also treats the normal conducting case. 

The thermal conductance is calculated on the level of 
the Kubo-Formula introduced in Ref. 83 for both normal 
metals and superconductors. In the conductivity bubble, 
the electrical charge in each current vertex is replaced by 
the half sum of adjacent fermionic Matsubara frequen¬ 
cies. Within the superconducting, interacting NLcrM ap¬ 
proach (see Sec. IIC 2), this amounts to introducing an 
energy dependent vector potential as a source field 84 


S a = 77 . f ti[(D,Q) 2 ]-2Z U f tr[(e + £A T v L afi )Q]. 
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The covariant derivative is 
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For the calculation of electrical conductivity, one replaces 
Z“ -*■ eJ“, where e is the electrical charge. Instead, for 
the calculation of thermal conductivity, we introduce 

(X“ 0 0 )l = 6 a ° a 5 a °P8 n - m>no i e -A^A. (B3) 


To obtain the thermal conductance including all quan¬ 
tum corrections stemming from scales L e ( 1,Lt c ) we 
perform RG up to the infrared cut-off scale and then 
evaluate the gradient term of the NLcrM in the saddle 
point approximation Q = A. We use 
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to determine the effective action for the source fields of 
the heat current 


S heat [A“J = -|E Al n Al_ n 
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For comparison, we present the same formula for the 
source fields of the electric current 


S e Utn\ = -jT j e 1 2 Al n Al_ 
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which can be evaluated to obtain the Mattis-Bardeen 
ac-conductivity formula in the superconductor includ¬ 
ing all quantum corrections stemming from scales L e 
(7, min {L u ,Lt c )}- When evaluating this expression in 
the limit A -» 0 we obtain 


c _ V 1 M acl 

be[ ~ E 2tt 2T Ai ’ n i ■" 


ot,n 


(B7) 


Now let us consider the superconducting state. First, 
we note that the sum in Eq. (B9) is formally UV di¬ 
vergent. This is a well known property of the Kubo ex¬ 
pression for the thermal conductance in superconductors. 
This divergence is cancelled by terms stemming from the 
time derivative of time ordering Heaviside functions (see 
e.g. discussion in Ref. 85) which is equivalent to disre¬ 
garding the large £fc contribution in the contour integra¬ 
tion. 

We define 
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and evaluate the sum S(ui) entering k = —5(w+zO)| w _ > o 


i.e. the normal state clc conductivity of cr = ge 2 /2n = 
ge 2 /h. 

We proceed with the evaluation of the thermal con¬ 
ductivity k entering the effective action of source fields 
as 
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and thus being for the superconductor (for simplicity we 
concentrate on ui n > 0) 
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where we analytically continued ui n ->■ ui + i0. We thus 
recovered the Wiedenrann-Franz law in the normal state 
(but including renormalization 81 ). 
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It follows that 
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Expanding in w we obtain Eq. (92) from the main text. 
On the Drude level for electrons with quadratic disper¬ 
sion we have g/2n = nr/m, and thus our result reduces 
to the formula of Ambegaokar and Griffin. 86 
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